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ABSTRACT 

We present an alternative interpretation of the observations of BR 1202—0725 at z = 4.695 and 
show that its properties are consistent with a relativistic jet, issuing from the quasar core, propagating 
into the inter-galactic medium and triggering star formation along its path. Prompted by this finding, 
we reviewed all the 2; > 3 objects detected in molecular line emission and found that the distribution 
of gas and dust in these sources is often spatially or kinematically offset from the host galaxy, and 
preferentially aligned along the radio axis. These observations suggest to us a scenario in which CO 
emission observed in high redshift galaxies is located where it forms: along the sites of star formation 
triggered initially by relativistic jets. 

Subject headings: galaxies: high redshift — galaxies: jets — stars: formation — early universe 



1. INTRODUCTION 

Ever since the first discoveries that high redshift 
quasars c an be brig ht sub-millimetre (submm) sources 
ijMcMaho n et alJT9 94') which shine because of dust emis- 
sion reprocessed from starlight, a large observational ef- 
fort has been undertaken to find corresponding reser- 
voirs of gas which provide both t he fuel and the site 
for this star forma tion to occur l)Omont et alJ Il99fit 
IGinlloteau etTTH TOTT. 1999). The gas reservoirs are as- 
sumed to be made up almost entirely of molecular hydro- 
gen (H2) with trace amounts of heavier molecules like 
Carbon Monoxide (^^C^^O; hereafter CO). But when 
and how did these molecules form? At high redshift 
it remains unclear how to collapse a first generation of 
metal-free stars and enrich the pristine primordial gas 
with metals and, eventually, molecules rapidly enough 
to approach solar metallici ty in quasars only 850 m illion 
years after the Big Bang ( Bromm fc Larso'iill2004(l . 
At z = 2, the co-moving space density of radio galax- 
ies with powerful jets was 1000 times higher than at 
the present epoch (^Willott et al. 200r). This hints to- 
ward the importance radio jets play in the early uni- 
verse. Between 2 < z < 3 there is a dramatic change in 
the optical morphology^ of the hosts of powerful radio 
galaxies, from massive ellipticals with old stellar popula- 
tions to irregular and clumpy distributions aligned along 
the radio axes fMcCarthv et alJ IT987t iChambers et al.l 
119871 Ivan Breueel et alJ Il998|l This alignment effect 
further indicates some connection between the forma- 
tion of the host galaxy and the coUimated jets from the 
AGN. The mere existence of powerful radio sources re- 
siding in such primeval galaxies hints toward a paradigm 
in which supcrmassive black holes are in place before 
their host fully forms. Although the current paradigm 
for galaxy formation is one in which black holes grow co- 
evally with host galaxies in hierarchical formation scenar- 
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ios ijWhite fc Reeslll978t iKauffmann et"anil993D . each 
of the above considerations could be consistent within a 
framework where jet-triggered star formation plays a key 
role in the formation of some of the first stars and galax- 
ies. In this Letter, we show that the z = 4.695 quasar 
BR 1202—0725 shows evidence for large scale jet-induced 
star formation; despite no direct detection of its radio jet. 
We then review all z > 3 molecular gas observations to 
date and show that often the gas and dust is aligned 
with the radio emission; in the case of bona fide radio 
galaxies this alignment is along the axes of radio jets. 
These results have led us to propose a scenario in which 
CO emission observed at high redshift is located where 
it formed: along the sites of star formation triggered ini- 
tially by relativistic jets. We assume a flat A-dominated 
cosmology with Hq — 71 km s^^ Mpc^^ and f^A — 0.73. 

2. BR 1202-0725 

At z = 4.695, BR 1202-0725 is an optically bright 

Mb = —28.5) quas ar with broad emission lines 

Storrie-Lombardi et al.l[T996i|l and the first high redshift 
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quasar confirmed to host massive re s ervoirs of molecu - 
lar gas and dust l|Omont et alJll996l: lOhta et al.lll996D 
Its 1.4 GHz luminosity (~ 10^^ W Hz"^: iCarilh et alJ 
l2002al hereafter C02) is consistent with either massive 
star formatio n (~ 2800 M© year~^; Equation 21 of 
ICondonllT99l or a Fan aroff-Riley type I radio galaxy 
ijFanaroff fc Riievl IT97|) . Luminous (~ lO^^ ergs s'^) 
X-ray emission from the quasar core provides clear evi- 
dence that BR 1202-0725 plays host to a powerful AGN 
(I. J. Klamer, D. M. Alexander fc C. Vigna l i, in p repa- 
ration). This is consistent with lYun et al.l l)2000D who 
measure a radio spectral index^ of —0.43 ± 0.14, which 
is flatter than the narrow range exp ected for pure star 
formation: {a) = -0.75, = 0.10 ()Condonl 119831 but 
consistent with a flat spectrum radio loud AGN. High 
resolution VLA observations, shown in Figure show 
two radio sources separated by about 24 kpc along a po- 

^ The spectral index, a (where Sv oc v") is measured from the 
integrated emission from BR 1202—0725 and its nearby compan- 
ion between 1.4 and 4.9 GHz. Future high resolution observations 
should be able to determine a separately for the two components. 
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sition angle (PA) of — 35±1.0°. The northern source, 
almost twice as luminous as the quasar at 1.4 GHz, coin- 
cides with a dusty, star forming companion (C02) which 
is undetected at optical wavelengths. If the radio emis- 
sion from the companion is dominated by star forma- 
tion, the corresponding star formation rate is a titanic 
~ 4400 Mq year~^; more than three times the rate im- 
plie d from the submm emission (Equation 26 of lCondonl 
llOg Z*). However, we find evidence to suggest that there 
is a relativistic jet emanating from BR 1202—0725 along 
a PA of about —35°. If correct, a portion of the radio 
emission could be due to synchrotron radiation from the 
relativistic electrons in the jet and the companion itself 
the result of jet-induced star formation; hereafter we refer 
to the quasar as 1202 A and the star forming companion 
as 1202B. This scenario would be somewhat analogous 
to the 850/zm SCUBA source LE 850.7 associated with 
j et-induced s tar formation from the radio lobe of R072 
pvison et al.l l2002'). 

Evidence for star formation occurring along a path from 
1202A to 1202B has been found in published optical and 
infra-red observations; these observations are depicted in 
Figure n Lyman— a emission (left panel; Figure from 
1202A is extended along a PA of — 28±5° with a velocity 
gradient shifting redward with increasing distance away 
from the quasar (Fontana ct al. 1998; Hu et al. 1996; 
hereafter H96). HST I-band observations (H96) reveal 
a faint, extended continuum source located about half 
way between 1202 A and 1202B along a position angle 
(PA) of -35±3° and 0'.'6 (4 kpc) from a Lyman-a emit- 
ter at z = 4.702. NICM0S3 K-band observations (right 
panel; Figure 13) of the same field (H96) shows that the 
quasar is clearly extended along a similar North- West 
direction. Two more K-band features are also detected, 
one located just North of 1202B along a PA of -36±3°, 
the other East of 1202B along a PA of -30±3°. The 
presence of such features, without corresponding I-band 
counterparts, is plausibly due to [Oii]A3727 line emission, 
a clear star formation indicator, which falls in the K-band 
window at this redshift. In addition, the CO line emis- 
sion is extended and the velocity profiles have structure. 
VLA imaging (C02) with 2.6" x 1.8" resolution reveals 
unresolved C0(2-l) emission at both 1202A and 1202B, 
whereas higher resolution observations (0.3" x 0.21") only 
detect 1202A confidently, suggesting 1202B has been re- 
solved. The CO line profile at 1202A is both stronger and 
narrower than at 1202B. Lastly, X-ray emission coinci- 
dent with 1202B has recently been detected (I. J. Klamer, 
D. M. Alexander & C. Vignali, in preparation) and is 
consistent with inverse Compton (IC) radiation against 
the cosmic microwave background (CMB). 
Although existing radio observations of 1202 A do not 
have the sensitivity to detect a jet, a scenario in which 
an inferred jet emanating from 1202 A is triggering star 
formation along its path toward 1202B, is a natural and 
consistent interpretation of all existing observations. The 
broad line profile of CO at 1202B would then refiect a 
bulk, turbulent velocity field associated with star for- 
mation from multiple cloud collapse and the gradient in 
the Lyman-a velocity would be consistent with ionisation 

along the axis of the jet. 

3. HIGH REDSHIFT CO OBSERVATIONS 

After finding evidence for large scale jet-induced star 
formation in BR 1202—0725, we collected the available 



data on all z > 3 CO emitters to investigate whether this 
model might be applicable elsewhere. Table 1 presents 
a brief summary of each source. The optical redshift 
is listed in Column 2 followed by the CO redshift in 
Column 3. Columns 4-7 summarise the radio, CO and 
submm information available for each source and the 
references used in constructing this table are listed in 
Column 8. Superscripts in Column 1 denote the source 
classification. Nine out of the twelve known z > 3 CO 
emitters show spatial or kinematic structure. Some char- 
acteristics are: 

• The CO and/or dust are spatially extended and aligned 
with the radio axis (e.g. B3 J2330+3927, 4C 60.07, 
6C 1909-^722). 

• The CO profile shows two kinematic and/or spatially 
distinct components and the host galaxy coincides with 
one (e.g. BRI 1335-0417, BR 1202-0725, 4C 60.07). 

• The CO velocity is offset with respect to the optically 
derived velocity (e.g. TN J0121-H1320, BRI 0952-0115, 
B3 J2330+3927, 4C 60.07, BRI 1335-0417). 

• The profile of any CO line detected at the position 
of a host galaxy is narrow compared with much broader 
CO emission offset from the host (e.g. BR 1202-0725, 
4C 60.07). 

Apart from the bona fide radio galaxies, there is still 
no direct evidence for jets in the other CO emitters. 
However, each of these properties is consistent within 
a framework in which the molecular gas and dust are as- 
sociated directly with star formation along the axis of an 
inferred relativistic jet. In this scenario, one would ex- 
pect that the CO profile from emission outside the host 
galaxy will be broad reflecting the bulk gaseous veloc- 
ity field from multiple cloud collapse, that the CO and 
Lyman-a will be tracing very different physical regions 
and so the derived redshifts for each will be different, 
and that the CO and dust will be extended and aligned 
with observable synchrotron radio emission. If jet-cloud 
interactions are responsible for CO production at high 
redshift, then many CO emitters should show evidence 
for jets. Out of the 12 2 > 3 CO emitters'^, six have 
evidence for jets aligned with either the CO or dust, 
five have radio luminosities above 10^^ W Hz~^ and are 
clearly AGNs, and a further four have radio luminosities 
above 10^^ W Hz~^ indicating either extreme starbursts 
or possible AGN. Figure |31 shows a histogram of relative 
position angles for the subset of CO emitters with re- 
solved radio and submm emission. Notwithstanding the 
small numbers, the clear excess centred around PA^ 0° is 
consistent with an alignment effect in z > 3 CO emitters. 
We note that in the absence of radio jets, an alignment 
effect is expected for star forming galaxies since the radio 
emission will be tracing the same physical regions as the 
gas and dust. However, all six sources considered in this 
histogram have AGN and at least three are radio galax- 
ies with detectable relativistic jets. In addition, there 
are many reasons why a jet may not be detected (or de- 
tectable) at radio wavelengths including (i) the active jet 
lifetime may be shorter than the period of continued star 
formation along the jet direction, (ii) radio observations 

^ Five of the sources are likely to be gravitationally lensed, with 
the potential for distortion and differential magnification caused by 
the foreground lens. Therefore, although we include these sources 
in our discussion, we note that their observations cannot be so 
easily interpreted. 
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made at rest frame centimetre wavelengths rather than 
(deci)metre wavelengths where lobes and hotspots are 
brighter, (iii) smaller magnetic fields or relativistic bulk 
flow with no relative motion between the electrons and 
the magnetic field, and (iv) rapid quenching of the radio 
emission due to IC losses against the CMB. Finally, since 
the CO emitters are heterogeneously selected (Table 1), 
they may not necessarily represent typical star-forming 
regions in the early universe. 

4. DISCUSSION 

4.1. The radio- optical alignment effect 

There is tantalising observational evidence which sup- 
ports jet-induced star formation as the natural explana- 
tion for the radio-optical alignment effect discussed in 
^ Well known examples exist in the local universe 
^n BreuEcl ct al. 1985; Rciku ba et al]l2002tl . but the 
evidence is most conspicuous in high redshift radio galax- 
ies like 4C 41.17 a.t z = 3.8 where unpolarised, rest- 
frame UV continuum is aligned along the radio axis 
l|Chambers et abllTflM IDev et aJjITflfll. If high red- 
shift star formation is first triggered on large scales due 
to the above process, then the first molecular gas and 
dust to form will be located in those regions where the 
first stars expelled their enriched material into the inter- 
stellar medium (ISM). Therefore, we predict that if CO 
is detected in 4C 41.17, it will be ahgned along the radio 
axis. 

4.2. In the beginning there were jets... 

Although the initial formation mechanisms of super- 
massive black holes remains largely unknown, the notion 
of seed black holes which form primordially and grow into 
a distribution of black hole mass es has been around fo r 
three decades (jCarr fc Hawkingll 974: Sil k fc ReesllOQ^ . 
The mass distribution would necessarily be governed, at 
least partially, by the density of the surrounding gas; 
the most massive black holes would then form in re- 
gions of the highest gas density and it will be these 
sites where we observe high redshift radio galaxies and 
radio-loud quasars. The highly relativistic, supersonic 
jets which power into the surrounding medium and slam 
into the existing overdensities can t rigger star formation 
along cocoon s surrounding the jets l)Bicknell et alJl2(¥)n : 
iFragile et al.l 12004) or could modify the stellar initial 
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mass function due to the effect of enhanced cosmic ray 
ionisation in the molecular cores (RDE; in preparation). 
Therefore, it will also be along these same preferential di- 
rections that the first heavy elements — including carbon 
and oxygen — will be produced as the stars end their lives 
and enrich their surroundings. This model provides the 
means of orchestrating star formation over tens of kpc 
on light crossing timescales. It is this process which has 
already been invoked to ex plain the ra dio-optical align- 
ment effect at high redshift (Rees 1989). Here we extend 
the idea to the formation of the first metals and molecules 
and predict, therefore, that a strong alignment exists 
generally between high-redshift molecular gas, dust and 
jets from supermassive black holes. 
But are the timescales consistent? Existing models for 
jet-induced star formation predict cloud collapse to oc- 
cur on timescales of order 10^ years. There will be a 
delay before CO formation takes place which, on the 
simplest level, will be the time it takes for the enrich- 
ment of the ISM. For 4C 41.17 (see gjl), the age of 
the stellar population aligned with the radio axis is a 
few X 10^ years. The lifetime of radio jets is typically a 
few X 10 years so multiple generations of stars can form 
along the axes of relativistic jets and consequently the 
CO to II2 abundance will increase to values approach- 
ing that seen in nearby galaxies. Once there are enough 
metals, molecules and dust available, then cooling and 
conventional star formation would take over, and be the 
dominant process in the local universe. 
Our hypothesis that AGN activity and jet formation pre- 
cede the assembly of their host galaxy has wide ramifi- 
cations. The continuum radiation from the jets could 
provide the fiux needed to obse rve Hi in absorption be- 
yond the epoch of reionization l|Carilli et al.ll200'2b|) . A 
jet-induced star formation scenario could also have pro- 
found implications for galaxy formation models. Should 
we be considering a paradigm in which the presence of 
such jets plays a role in deciding the morphology of the 
host galaxy? 

We thank C. Carilli and E. Hu for the images of 
BR 1202—0725, D. Alexander for valuable discussions 
and the anonymous refereee for useful comments which 
improved this Letter. EMS and RWH acknowledge grant 
support from the Australian Research Council. 
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TABLE 1 

Radio, CO and submm properties of all published z > 3 CO emitters. 



Source Zopt ZCO Radio PA Radio Power CO PA ^ Dust PA ^ References 









O 


W Hz-i 


O 


O 




SDSS J1148+5251'' 


6.41 


6.42 




1 X 10^^*^ 






1,2,3,4 


BR 1202-0725'^ 


4.694 


4.695 


-35 


3 X 10^5 


-35 


-40 


5,6,7 


BRI 0952-0115='^ 


4.426 


4.434 




4 X lO^Sg 






5,8,9 


BRI 1335-0417'^ 


4.396 


4.407 




3 X 10^^ 


+13 


+18 


5,8,9,7 


PSS J2322+1944<^ 


4.111 


4.119 





4 X 102**s 


+5 




10,11,12 


APM 08279+5255<^ 


3.87 


3.91 


+23 


1 X lO^Sg 


+20 




13,14,15 


4C 60.07° 


3.788 


3.791 


-73 


4 X 10^* 


-75:'^ 


-70 


16,17,18,19 


6C 1909+722<= 


3.536 


3.532 


+15 


5 X 1028 




+18 


20,17,18 


TN J0121+1320° 


3.516 


3.520 




1 X 1028 






20,21 


SMM 09431+4700^^ 


3.35 


3.35 




4 X 10^4 fg 






22,23 


B0751+2716'*'= 


3.200 


3.200 




1 X lO^Sg 






24,25 


B3 J2330+3927° 


3.087 


3.094 


+33 


9 X 10^7 




+50 


26,18 



References. — (l)l\ViUott ot aTllaOOa ( 2) IWalter et al.ll20ci^ (3) |Bert oldi ct al. 20q3 | (41 [C arilli et al.ll2004l (5) IStorrie-Lombardi et _al.lll996i (6) 
gmonyit ai. 1996 (7) Carilli ct al. 2002a (8) Guillotcau ct al. 1997 (9) Guillotcau ot al. ' 1999 (10) ^CMillLfit_aI.'"2001 (11) Carilli ct al. 20 03(12) 
Cox et a l. 200 2 (13) Ib ata ct al. 1999 (14) Papadopoulos ct al. 2001 (15) Downcs ct al. 1999 (16) C ha mbers et al. 1996 (17) PapadoDouloTe^U 
200a (18) Stcvens^^^lJ|200^19) |Greve_e^^lJ|200|L(20) |De_Breuck_e£ay|2001|_[21) IDe Breuck et al.ll2003al (22) fCowie et al.ll2002i (23) jNe^^Q 
200 3 (24) B arramisTral1l200a (25T lLehar7r^lT99^ ( 26T lDe^reuckeraUI2003bl 



Note. — ^Radio luminosity measured at 1.4GHz in the rest frame of the source. PAs are measured from published images given in the reference 
list. Errors in PA are of order 5°. ^Optically selected quasar with submm emission. Gravitationally lensed with submm emission. ^Radio galaxy 
with submm emission. ^Assuming a radio spectral index of —0.7. ^Not corrected for possible gravitational magnification. ^Based on the broad 
CO(l-O) component of lGreve et al.l I2004H . 
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Fig. 1. — 1.4 GHz radio continuum contours of BR 1202—0725 with pluses marking the C0(2-l) positions (C02). Contour 
levels are a geometric progression of starting at ±42 A ^Jy beam~^. Striped, bricked and filled circles depict K-band, I-band 
and Lyman-a features respectively. North is up and East is to the left. A bar in the top left corner shows the 2" (13.2 kpc) 
scale. 




Fig. 2. — Multi- wavelength observations of the quasar BR 1202—0725. Left: Lyman-a greyscale with CO(2-l) contours overlaid 
(C02). Right: NICMOS3 K-band emission from H96. The open circles correspond to the K-band features indicated in FigureQ 
North is up and East is to the left. A bar in the top right corner of each panel show the 2" (13.2 kpc) scale. 
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Fig. 3. — Stacked histogram of relative differences between the radio-CO (shaded) and radio-dust (unshaded) PAs for the 
six sources in Table 1 with sufficient PA information. BR 1202—0725 and 4C 60.07 have both radio-CO and radio-dust PA 
information and are counted once in the shaded region. In the absence of an alignment effect, the relative PAs will be randomly 
distributed into nine bins over 0° <PA< 90° and the probability that five of the six objects fall in the 0° — 10° bin is 1.5 x 10~^. 



